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To analyze the physical properties arising from indirect magnetic exchange between several mag- 
netic adatoms and between complex magnetic nanostructures on metallic surfaces, the real-space 
extension of dynamical mean-field theory (R-DMFT) appears attractive as it can be applied to 
systems of almost arbitrary geometry and complexity. While R-DMFT describes the Kondo ef- 
fect of a single adatom exactly, indirect magnetic (RKKY) exchange is taken into account on an 
approximate level only. Here, we consider a simplified model system consisting of two magnetic 
Hubbard sites ( "adatoms" ) hybridizing with a non- interacting tight-binding chain ( "substrate sur- 
face"). This two-impurity Anderson model incorporates the competition between the Kondo effect 
and indirect exchange but is amenable to an exact numerical solution via the density-matrix renor- 
malization group (DMRG). The particle- hole symmetric model at half- filling and zero temperature 
is used to benchmark R-DMFT results for the magnetic coupling between the two adatoms and 
for the magnetic properties induced in the substrate. In particular, the dependence of the local 
adatom and the nonlocal adatom-adatom static susceptibilities as well as the magnetic response of 
the substrate on the distance between the adatoms and on the strength of their coupling with the 
substrate is studied. We find both, excellent agreement with the DMRG data even on subtle details 
of the competition between RKKY exchange and the Kondo effect but also complete failure of the 
R-DMFT, depending on the parameter regime considered. R-DMFT calculations are performed 
using the Lanczos method as impurity solver. With the real-space extension of the two-site DMFT, 
we also benchmark a simplified R-DMFT variant. 

PACS numbers: 71.10.Fd, 71.27.+a, 75.20.Hr, 75.75.-c 



I. INTRODUCTION 

The rapidly improving experimental techniques to 
probe magnetic adatoms on non-magnetic surfaces al- 
low for direct studies of fundamental magnetic exchange 
mechanisms on an atomic scale. Besides access to the 
structural and the electronic properties of such adatoms 
and of the underlying substrate for a given system, the 
construction of tailored magnetic model systems repre- 
sents an exciting perspective^— Magnetic structures of 
nanometer size provide extremely small systems suitable 
to store and to transport information and may realize 
efficient nano spintronics devices £ 

The competition between an indirect magnetic ex- 
change of the adatoms via the substrate electrons on 
the one hand and the screening of the adatom mag- 
netic moment by the conduction-band electrons of the 
substrate on the other represents a prominent exam- 
ple for a physical problem becoming accessible to new 
real-space techniques. The scanning tunneling micro- 
scope (STM)& has been used to investigate the Kondo 
physics! of single magnetic adatoms^^ and the mag- 
netic properties of the individual magnetic island o 10 i n 
on non-magnetic substrates. Using STM, it is pos- 
sible to investigate the direct magnetic interaction of 
atom pairs ' 12 ' 13 Indirect magnetic exchange, i.e. the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction^ 
between two adatoms, was detected through the Kondo 
effect.— A direct real-space study of the RKKY coupling, 
however, comes in reach with spin-polarized scanning- 
tunneling spectroscopy onlyJ£~— 

The most simple model which captures this compcti- 
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FIG. 1: (Color online) Schematic picture of real-space dy- 
namical mean-field theory (R-DMFT) for the two-impurity 
Anderson model (TIAM). The system is given by two "mag- 
netic" sites with strong Hubbard interaction U (orange) at a 
distance d coupled via a hybridization term of strength V to 
a one-dimensional "substrate" consisting of L non-interacting 
sites (blue) with nearest-neighbor hopping t = 1. In the R- 
DMFT, the TIAM is self-consistently mapped onto two single- 
impurity Anderson models which are solved independently 
by means of exact diagonalization (ED) to get the local self- 
energies. These are used to set up the TIAM Dyson equation 
the solution of which gives the local Green's functions which 
define via the R-DMFT self-consistency conditions the pa- 
rameters of the impurity models (see text for details). 



tion is displayed schematically in Fig. [T] Here the elec- 
tronic and magnetic properties of a magnetic adatom 
are modeled by a single non-degenerate orbital. A lo- 
cal magnetic moment is formed by a strong local Hub- 
bard interaction. The adatom orbital hybridizes with 
a valence orbital of the nearest-neighboring substrate 
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atom. The substrate electronic structure itself is mod- 
eled by a tight-binding valence band resulting from non- 
degenerate and uncorrelated orbitals on a bipartite lat- 
tice with nearest-neighbor hopping. Considering two 
adatoms yields a variant of the two-impurity Anderson 
model (TIAM)i£ in a surface geometry. The main goal of 
our study is to benchmark the real-space variant of the 
dynamical mean-field theor y 21 ' 22 that can be employed 
for theoretical studies of the electronic and magnetic 
properties of a single, of two and or of more magnetic 
adatoms in different geometries on metal surfaces. For 
this purpose the TIAM represents a fundamental starting 
point. In the Kondo limit of the TIAM, charge fluctua- 
tions on the adatom site are largely suppressed, and the 
adatom spin couples antiferromagnetically to the lo- 
cal spin at the nearest-neighboring substrate site S c via a 
spin-spin coupling —JS^S C given by the local exchange 
Joe -V 2 /U < 0. 

The interplay between the Kondo effect and the 
RKKY interaction has extensively been studied in the 
Kondo limit of the TIAM or in the two-impurity 
Kondo model by different analytical as well as numer- 
ical techniques ! 19 ! 23 " — The physics is governed by two 
energy scales, the nonlocal indirect magnetic interac- 
tion Jrkky oc J 2 and the Kondo temperature Tk oc 
exp(—l/\J\) below which the magnetic moment of the 
adatom (impurity) is screened locally. In the Kondo 
regime for Tk ^> | Jrkky |, the conventional picture is 
that the local magnetic moments at the two impurities 
are individually screened by forming local singlet states 
with two Kondo clouds of itinerant electron spins from 
the substrate (conduction band). For large | Jrkky |, 
on the other hand, and in the antiferromagnetic case 
Jrkky < 0, the two adatom spins form a nonlocal sin- 
glet state and there is no Kondo effect. If Jrkky > is 
ferromagnetic and large as compared to Tk, a nonlocal 
spin-triplet state is formed. This may subsequently be 
Kondo screened. In the generic case and as a function 
of J there is no quantum phase transition but a smooth 
crossover from the RKKY regime at weak J to the Kondo 
regime at strong J. For a dense system, i.e. the Kondo 
or Anderson lattice model, a static mean-field approach 
would sharpen this to a phase transition^ 

Dynamical mean- field theory is a comprehensive, ther- 
modynamically consistent and non-perturbative approx- 
imation for correlated lattice-fermion models^ DMFT 
treats the Kondo effect exactly. On the other hand, 
one has to tolerate an approximate treatment of the ef- 
fects of the RKKY interaction. It is important to note 
that there is no approximation of RKKY coupling it- 
self: Integrating out the non-interacting substrate de- 
grees of freedom, the effective second-order-in- J RKKY 
coupling, Jrkky, ij = J 2 X°/ U (w = 0), is given in terms 
of the nonlocal static susceptibility of the substrate. It 
was pointed out by Peters and Pruschke^ that Jrkky, ij 
is still finite but reduces to an interaction between near- 
est neighbors for the case of a lattice in infinite spatial 
dimensions where the DMFT becomes exact. For finite 



dimensions, it is a long-ranged and oscillating function 
of the distance d = \i — j\. What is neglected in fact 
for a finite-dimensional lattice, is the feedback of non- 
local, e.g. magnetic correlations, which result from the 
nonlocal RKKY coupling, on the self-energy and thus on 
the one-particle Green's function. This is a rather subtle 
approximation the quality of which can be estimated by 
concrete numerical calculations only. 

The same argumentation holds for the real-space 
DMFT (R-DMFT) 3 R-DMFT generalizes the standard 
DMFT to systems with missing or reduced translational 
symmetry by self-consistently mapping the original (lat- 
tice) model to a set of single-impurity Anderson models 
(SIAM), one for each of the geometrically or electroni- 
cally inequivalent sites. Even for the TIAM (see Fig.[T|), 
this real-space generalization is necessary if one wants to 
apply DMFT in order to test the local approximation for 
the self-energy. 

Previous applications of the R-DMFT concentrated 
on the Mott metal-insulator transition at surfaces and 
in thin filmS ) 20 i 46 on surface effects in correlated Fermi 
liquids," on multilaycred nanostructurcs, heterostruc- 
tures and interfaces^ on disordered systems r 9 - as well 
as on ultracold atomic gases in optical lattices with har- 
monic confinement^ It has not been employed, however, 
to study the effects of the indirect magnetic exchange. 

The main purpose of the present study is to apply 
the R-DMFT to the particle-hole symmetric TIAM at 
half-filling and zero temperature and to study the mag- 
netic response, i.e. different static magnetic susceptibili- 
ties, by applying a weak local magnetic field to one of the 
adatoms. Calculations are performed as a function of the 
distance between the adatoms and as a function of the 
hybridization strength V to cross over from the Kondo to 
the RKKY regime. To test the reliability of the dynam- 
ical mean-field approach, the substrate electronic struc- 
ture is modeled as a one-dimensional tight-binding chain 
(see Fig. The resulting essentially one-dimensional 
model is accessible to the density-matrix renormalization 
group (DMRG)i£i~— Extensive comparison with numer- 
ically exact DMRG results obtained from an implemen- 
tation based on matrix-product states^ and operators 
along the lines described in Ref. [55] helps to benchmark 
the mean-field approach. 

Our intention is that, by comparing with DMRG, the 
strengths but also the mean-field artifacts of R-DMFT 
become more transparent. A failure of R-DMFT for the 
weak-coupling limit, where non-local correlations due to 
the RKKY coupling are strong, can be expected from the 
very beginning. However, there are several interesting 
questions left, e.g.: Where precisely are the limits of the 
mean-field approach? How does a failure of the approach 
manifest itself in the observables? Which physical effects 
are accessible to a description on the R-DMFT level? 
To what extent can the physics be reproduced quantita- 
tively in the strong-coupling limit? Such benchmarking 
of the R-DMFT, at the level of the two-impurity Ander- 
son model, will be important for future studies of similar 
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systems in higher spatial dimensions, with more corre- 
lated adatoms forming more complex geometries such as 
chains or clusters etc. By choosing the one-dimensional 
two-impurity Anderson model at half-filling, the above- 
mentioned questions are tackled in a situation that is 
very unfavorable for R-DMFT. The benchmark will thus 
serve as a "lower limit" for the applicability of R-DMFT 
for future applications. 

Our interest in the R-DMFT approach to study mag- 
netic nanostructures on surfaces results from its ex- 
tremely large flexibility. Opposed to DMRG, for ex- 
ample, the R-DMFT is able to investigate inhomoge- 
neous systems in arbitrary geometries in higher dimen- 
sions. While this is actually characteristic for any mean- 
field approach, the R-DMFT is distinguished by the fact 
that it is non-pcrturbative and thermodynamically con- 
sistent. To account for the effects of short-range correla- 
tions, the theory can be improved by certain cluster ex- 
tensions, such as cellular DMFT.— This is, in principle, 
also conceivable for complicated inhomogeneous geome- 
tries but requires further methodical advances as there is 
no straightforward tiling of the lattice in most cases. 

The paper is organized as the follows: The next sec- 
tions introduces the model, notations and quantities of 
interest. Sec. Mil and Sec. IIVI briefly describe our real- 
space DMFT and our DMRG approach to the problem, 
respectively. Results of both approaches are presented, 
compared and discussed in detail in Sec. [V] Finally, Sec. 
IVTl concludes the paper. 



II. MODEL AND BASIC THEORY 

The Hamiltonian of the two-impurity Anderson 
models displayed in Fig. [T]is given by: 
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Here f\ a and cj a create an electron with spin projec- 
tion a =t, 4- a t the adatom sites a = 1, 2 or at the sub- 
strate sites i = 1,2, ... ,L, respectively. n^ a a = fa, a .f a ,cr 
and n c ia = cj rT c i a denote the corresponding occupation- 
number operators. The spin-summed occupation at one 
of the adatom sites and at one of the substrate sites are 



given by n£ = W a ^ + n' ^ and n-f = n" i1; + n\^, respec- 
tively. The hopping amplitude t between neighboring 
substrate lattice sites is used to fix the energy unit, i.e. 
t = 1. V is the hybridization between an adatom site a 
and the nearest-neighbored substrate lattice site which is 
denoted by i a . U and e are the on-site Hubbard inter- 
action and the local on-site energy for the adatom sites, 
/i is chemical potential. In all our calculations we con- 
sider the particle-hole symmetric case with fi = and 



e = —U/2 where the system is half- filled, i.e. where the 
average occupation numbers in thermal equilibrium are 
given by (n£) = 1 and (nf) = 1 for both a and all i. 

The magnetic properties of the system are best char- 
acterized by site-dependent local and nonlocal suscepti- 
bilities. We consider the adatom-adatom susceptibilities, 
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i.e. the local adatom susceptibilities Xaa for a = 1,2 
and the inter-adatom susceptibility X12 = X2i- These 
provide information on the local adatom magnetic mo- 
ment and, most important, on the indirect magnetic 
coupling. Further, we are interested in the linear mag- 
netic response of the substrate which is accessible via the 
adatom-substrate susceptibilities 
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Here m£ = (S^J and m? = (SfJ, with S^ z = |(n^ f - 

n a i) an d ^tz — ^( n i t ~ n i |)' are magnetic moments on 
the adatom site a and on the substrate lattice site i re- 
spectively. Furthermore, the imaginary-time dependence 
of an operator A is given by A(t) = e Hr Ae~' HT In our cal- 
culations the susceptibilities, Eqs. (J2|) and ©, are com- 
puted as a numerical derivative with respect to a local 
magnetic field of strength hp coupling as % — > 'H — hpSp^ 
to the Hamiltonian. Calculations are done for zero tem- 
perature in the present study. Nevertheless, the formal- 
ism is set up for arbitrary finite T below. 

The main task is to compute, for a finite but weak 
field hp, the spin-dependent average occupation numbers. 
These can be obtained via 



« CT ) = ^ + 2TEReG^y;<) 
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n>0 



from the local single-electron adatom Green's function 
G [ Z P A^n) = ({f a y,fl, a )) Un and the local substrate 
Green's function G^(iuj n ) — ((c, l a ; c\ )) u „ given at the 
fermionic Matsubara frequencies u> n = {2n + l)nT. The 
local Green's functions are the diagonal elements of the 
Green's function matrix G a (iuj n ). The latter can be ob- 
tained from the real-space Dyson equation: 



G^^Un) = (iuj n +fi)I — e cr -T- S CT (io; n ) 



(6) 



where / is the unity matrix, e a the diagonal local en- 
ergy matrix, and / is the hopping matrix. e a also in- 
cludes the field term and is thus possibly spin-dependent. 
T not only includes the hopping t between substrate 
sites but also hopping V between the substrate and the 
adatom sites. For a system with L substrate sites and 
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two adatoms, the matrix dimension is L + 2 for each spin 
direction a. 

As there is a local Hubbard interaction on the adatom 
sites only, the self-energy S ff (i«„) is a L + 2-dimensional 
matrix with non-zero elements £ a( g j(T (iw n ) in the 2x2 
adatom-sites block only. Hence, Eq. ^ can be written 
as: 
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(7) 



where Ca,<r = Ca,tr(iw n ) = zu; n + fj, - e a - £ QQ;(T (iu; n ), 
and where V is the 2 x L hybridization matrix including 
hopping between adatoms and substrate only. Its non- 
zero elements are given by V\a x = Vi,^ = V. Further, G° 
is the non-interacting substrate Green's function matrix. 
In case of periodic boundary conditions, its elements are 



L ^ iu n + (i- s{k m ) 

ra—0 v ' 



(8) 



where k m = 2irm/L with m = 0, 1, . . . , L — 1, while for 
open boundary conditions 



G 



sin(fc m z)sin(fc m j) 
+ 1 ^ ioj n +fi- e(k m ) 



(9) 



where k m = 7rm/(L + 1) with m = 1,2, . . . ,L. In both 
cases e(k) = — 2icos(fc) is non-interacting dispersion. 

Generally, to calculate, for a given self-energy, the local 
Green's functions, one has to numerically invert the ma- 
trix given by Eq. In our case for the TIAM, however, 
it is possible to "integrate out" the substrate degrees of 
freedom and to find analytical expressions which substan- 



tially reduce the numerical effort: Using the identity 
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with F = [A — BD 1 C) which is valid for arbitrary 



quadratic matrices A and D and arbitrary rectangular 
matrices C and D. we find 
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where 




is the 2x2 adatom Green's function matrix and 

Q S ub _ q0 _|_ qOv^^VG (13) 
is the L x L substrate Green's function matrix. The 
remaining task thus consists in the inversion of a 2 x 2 
matrix 



G imp (iuj ) = [ < ' 1,CT _ ^ n ~^12,<r — ^12 
I -S 2 i - A 2 i I<7 Ca.o-^n) - A 22 



where 
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(15) 



is the hybridization function. This is readily done: 



/-(imp / ■ \ (2,a(i0J n ) ~ A 22 (lW„) 
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r imp r ) = Cl,a{iu„) ~ Au(iu n ) 

^imp/. s _ ^yimp /• n_ Ai 2 (fcj„) + Si 2 , cr (lW n ) 

(Ci.<r( lw «) - An(za; n ))(C2, -(«a;„) - A 22 (iw„)) - (A 12 (ia;„) + £i2,<x(w<>nj) 



This provides us with the local adatom Green's functions in particular and, using Eq. (|13[) . with the local Green's 
functions for each substrate site via: 

G8£(«u„) = G&(iw n ) + £ G° iia {iu n )VG™» c {iu n )VG° ifj Sun) ■ (19) 



III. REAL-SPACE DYNAMICAL MEAN-FIELD 
THEORY 

To complete the theory, we need the self-energy ma- 
trix E Q ^ j0 -(iw n ). This requires an approximation. Within 



real-space DMFT^ the self-energy is obtained by con- 
sidering weak-coupling perturbation theory in U to all 
orders and by summing all local diagrams in the skeleton- 
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diagram expansion of the self-energy, S = £[G]. This im- 
plies that the resulting self-energy is local: £ Q( 3 jCr (iu; n ) = 
S a ^a,i7{iijJn) but possibly site-dependent. For correlated 
lattice models with full translational symmetries, the ap- 
proach reduces to the conventional DMFT— As in the 
conventional DMFT, the local diagrams are not summed 
explicitly, the problem is rather reformulated by intro- 
ducing a self-consistent mapping onto an effective single- 
impurity problem. Here, however, the self-consistent cy- 
cle is more complicated since a lattice model with M ge- 
ometrically or electronically incquivalent sites has to be 
self-consistently mapped onto a set of M effective single- 
impurity models. In our case we have to consider at most 
M = 2 single- impurity Anderson models (see Fig. [T|): 

We start with a guess for the local self-energies 
Sa,ff(jwn), i-e. for Ca,o-(*w n ). This is used in the Dyson 
equation of the lattice model to compute the Green's 
function matrix, and in particular the local elements of 
the Green's function matrix at the correlated sites. In 
our case, we can profit from Eqs. (flg)| and (fT7|) to get the 
local adatom Green's functions G™ p (T (icj„) directly. The 
R-DMFT self-consistency conditions, 

7^ — 7 r = — : h S Q a{ioJ n ) , (20) 

then provide us with the Weiss Green's functions 
£7° a (iuj n ) for a = 1, ...,M, i.e. with the non-interacting 
Green's functions of the M effective impurity models. 
These can be written as o-i^n) — * w n + ^ — £ — 
A a (iu) n ). The one-particle parameters of each effective 
SIAM, the one-particle energies of the bath sites as well 
as the corresponding hybridization strengths, are found 
from the poles and the residues of the corresponding hy- 
bridization function A a (iuj n ) (which should not be mixed 
up with A a p{iuj n ), see Eq. ([15)) ). Once the effective im- 
purity models are fixed, the crucial step consists in the so- 
lution of the models which can be done independently for 
any a = 1, M. This yields the self-energies Yi a ^ a (iu) n ) 
and thus closes the self-consistency cycle (Fig. [1} . This 
procedure is iterated until converged self-energies are ob- 
tained. 

As an impurity solver to get the self-energy E Q , jCr (iu; rt ) 
of the a-th SIAM we use the exact-diagonalization (ED) 
mcthod i 57 ' 58 Here, a finite small number n s — 1 of aux- 
iliary bath degrees of freedom in the effective SIAM is 
considered. We use full diagonalization with n s = 6 and 
the Lanczos metho d 59 ' 60 with n s = 8 and n s = 10. Ex- 
ploiting the fact that the total particle number and the 
z-component of the total spin are conserved quantities, 
the diagonalization can be done in smaller invariant sub- 
spaces of the full Hilbert space. All calculations have 
been done, if not stated differently, with n s = 10. For a 
given Weiss Green's function, the one-particle bath pa- 
rameters of the SIAM are found by a minimization pro- 
cedure on the imaginary-frequency axis as described in 
Ref. using high-frequency cutoff of the order of U and 
low-frequency cutoff specified by the fictitious tempera- 
ture T/t — 0.001. With the latter we can formally work 



in the finite-temperature Matsubara framework as out- 
lined above. On the other and, the value of the fictitious 
temperature chosen is clearly lower than the smallest en- 
ergy scale that can be accessed by means of the ED solver 
for n s = 10. We have regularly checked that the results 
do not significantly depend on the cutoff and on n s . 

The computational effort of the R-DMFT scheme 
roughly scales linearly with the number of impurity mod- 
els, i.e. with the number of inequivalent sites in the orig- 
inal system. While here we focus on the M = 2 case 
for benchmarking purposes, future applications are in- 
tended that address systems with up to 0(100) inequiv- 
alent magnetic atoms. We expect that those applications 
can still be performed conveniently using ED as a solver. 
It might nevertheless be interesting to have a scheme 
at hand that is considerably faster. Here, the two-site 
DMFT— represents an alternative. The main idea is the 
use a single bath degree of freedom only, as in the so- 
called linearized DMFT-S for the Mott-Hubbard tran- 
sition. The resulting effective two-site impurity model 
is readily solved. On the other hand, the DMFT self- 
consistency condition can obviously no longer be satis- 
fied exactly or to a high level of accuracy as in the ED 
approach with, say, n s = 10 sites. It has therefore been 
suggested 6 ^ to take into account the leading orders in sys- 
tematic expansions of the self-consistency condition for 
high and for low frequencies only. This results in a simpli- 
fied but extremely fast approach which is suitable to get a 
quick overview of magnetic phase diagrams, for example. 
The real-space extension of two-site DMFT to systems 
with reduced translational symmetries is straightforward 
and will be discussed in detail elsewhere. 



IV. DENSITY-MATRIX RENORMALIZATION 

The two-impurity Anderson model in an essentially 
one-dimensional geometry (see Fig. [IJ is amenable to a 
numerically exact solution by using the density-matrix 
rcnormalization group—"— Therefore, DMRG calcula- 
tions can be used to benchmark the quality of the mag- 
netic susceptibilities obtained from the R-DMFT ap- 
proach. The calculation of the ground state and of 
ground-state expectation values for the TIAM is a stan- 
dard problem within DMRG. Here we have been follow- 
ing Ref. [55] and have implemented a code which is based 
on the variation of matrix-product states (MPS)— The 
main idea is to optimize a test wave function \ip) of the 
form 

|V;> = A^...A^)|nx)...K), (21) 

ni,...,rj£ 

where is a local basis at the site q of a one- 

dimensional chain with L sites in total. The elements 
■^ilii °f the matrices A'™ 5 ' are considered as varia- 
tional parameters which are locally and iteratively opti- 
mized during a sweep through the chain by exploiting the 
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Ritz variational principle. Several sweeps are necessary 
to obtain a converged ground state. 

In practice, the local optimization can be reformulated 
as a generalized eigenvalue problem which is simplified to 
an ordinary one by exploiting a local gauge invariance 
of \ip) to properly (left- and right-) orthogonalize the 
A-matrices. The eigenvalue problem is then efficiently 
solved by means of the Davidson method^ We profit 
from the so-called wave-function transformation^ to re- 
duce the number of iterations necessary for convergence 
of the Davidson algorithm and exploit the two U(l) sym- 
metries of the Hamiltonian corresponding to conservation 
of the total particle number and the z-component of the 
total spin. 

Observables and the Hamiltonian in particular are rep- 
resented as matrix-product operators^ Besides an ele- 
gant and flexible coding this allows to easily consider 
different implementations of the Hamiltonian. For the 
present case of the TIAM there are two possibilities to 
treat the adatoms suggesting themselves: (i) An adatom 
orbital a and the substrate orbital i a "below" a are 
treated as a single "site" q in the DMRG context. The 
disadvantage is that therewith the local Hilbert-space di- 
mension at q is enlarged, (ii) The adatom orbitals a are 
treated as separate sites, i.e. a chain of length L + 2 is 
formed. This leaves the local Hilbert-space dimension 
constant but introduces next-nearest-neighbor hopping 
terms. We have tested both variants and found the differ- 
ences in computational costs and accuracy to be marginal 
only. Routinely, variant (ii) is employed. 

An important aspect is to prevent the sweep algorithm 
from getting stuck in a local energy minimum. This 
can be circumvented by implementing a mixed single-site 
approac h 55 ! 65 to introduce fluctuations in the reduced 
density matrix. The additional coupling to a larger set of 
states considerably improves the convergence properties 
when optimizing It furthermore also allows to dy- 
namically adapt the dimensions of the A-matrices. Con- 
verged results for typical situations with long-range spin- 
spin correlations in a TIAM with about L = 50 sites are 
obtained with matrix dimensions of the order of m = 400 
in the largest invariant blocks of the ^-matrices. A re- 
liable error measure is the variance r = (ipKH — E) 2 \ip) 
which is easily accessible within an MPS-based imple- 
mentation. We have checked that typically the standard 
deviation -Jr < 1 0~ 4 . 



V. RESULTS 

R-DMFT and DMRG calculations have been per- 
formed for the TIAM at half-filling and zero tempera- 
ture. We consider systems with an even number L of 
substrate sites and two adatoms at positions symmetric 
to the chain center at a distance d — \ii — i%\ as displayed 
in Fig. [TJ As this implies an even number of electrons, 
there is no Kramers degeneracy of the ground state. All 
calculations are done using open boundary conditions. 



Fig. [3] shows the magnetic susceptibility xTp 3 f° r a 
TIAM with L = 50 as defined in Eq. ©. Because of the 
mirror symmetry, it is sufficient to discuss e.g. the case 
(3 = 1, i.e. the left adatom. For the calculations we apply 
a weak local magnetic field with a strength hp=\ ~ 10~ 5 - 
10~ 2 at the left adatom and look for the response at 
substrate site i. Regularly, calculations for different hi 
are performed to ensure that the field strength is in the 
linear-response regime. 

Let us first concentrate on distances d = 4n+l with in- 
teger n. Here the RKKY coupling between the magnetic 
adatoms is antiferromagnetic. Other distances d includ- 
ing those with ferromagnetic coupling will be discussed 
in Sec. IVEl 



A. Magnetic response of the substrate 

We start the discussion with d = 49 (lowest panel 
on the right in Fig. [2]). This is the case where the 
two adatoms are located at the edges of the substrate 
chain. The blue lines refer to our R-DMFT calculations 
which have been done with n s — 10 local degrees of 
freedom in the effective impurity model. Directly "be- 
low" the first adatom at i\ the response is antiferromag- 
netic, i.e. Xii.i < 0. This simply reflects the antiferro- 
magnetic Kondo coupling J. The calculations have been 
done for U = 8 and V 2 = 2 where the nearest-neighbor 
hopping in the substrate t = 1 is used to set the en- 
ergy scale. This results in a negative, i.e. antiferromag- 
netic, local exchange interaction of intermediate strength 
J = —8V 2 /U = —2. This is clearly beyond the weak- 
coupling limit J — > but still charge fluctuations are 
largely suppressed: We find an average double occupancy 
of (711-f-ni^) = 0.072 at the adatom site, and the adatom 
local magnetic moment (Sf) = 3(1 — 2(ni^nii)) / 4 = 0.64 
is much closer to the localized-spin value 3/4 than to the 
free fermion value 3/8. 

As a function of the distance \i — i±\ to the first im- 
purity, the response is oscillatory corresponding to the 
2kp = 7r nesting wave vector. Its absolute value is max- 
imal at ii, decreases with increasing i and almost satu- 
rates until there is a slight upturn for i — > 12 = 50, i.e. at 
the position of the second adatom. Consistent with the 
2kp oscillation, xTi* 1S positive at i = which implies, 
due to the antiferromagnetic local coupling J < 0, that 
there is an antiferromagnetic (RKKY) alignment of the 
two adatom moments. 

The corresponding DMRG results are also shown in 
Fig.[5]for comparison (red lines). For the distance d = 49, 
however, there is actually no difference to the R-DMFT 
results visible on the scale of the figure. As R-DMFT 
accounts for the single-impurity Kondo effect exactly, 
this perfect agreement would be plausible if a picture of 
two independent Kondo effects applied. Strictly speak- 
ing, however, this cannot be the case: There is a finite 
nonlocal adatom-adatom susceptibility, even in this long- 
distance limit (see also Fig. [5] and corresponding discus- 
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FIG. 2: (Color online) Static magnetic susceptibility Xi" b a t Hubbard interaction U — 8 and hybridization strength V = 
for a system with L — 50 substrate sites and two adatoms at positions symmetric to the chain center and different distances d 
as indicated. Xii b gives the linear response of the substrate at site i to a local magnetic field at the first (left) adatom. Energy 
scale: nearest-neighbor hopping in the substrate t = 1. Blue lines with squares: results as obtained from real-space DMFT 
using exact diagonalization with n s = 1 as a solver. Red dashed lines with circles: numerically exact solution as obtained from 
DMRG calculations. For comparison DMRG results for a system with L — 49 sites and a single adatom are shown (dotted 
green line with triangles) . This corresponds to switching off the hybridization V between the substrate and the second (right) 
adatom. 



sion below) which in principle has a non- vanishing feed- 
back on the self-energy and generates nonlocal elements 
of the self-energy in particular. R-DMFT is thus approx- 
imate. On the other hand, we can conclude that this 
feedback of the nonlocal susceptibility is apparently neg- 
ligibly small and R-DMFT almost exact for the present 
situation. 

Upon decreasing the distance between the adatoms, 
this picture should change gradually. However, apart 
from the extreme case d = 1, deviations of the R-DMFT 
from the DMRG results arc extremely small, and the 
agreement between R-DMFT and DMRG remains excel- 
lent. On the other hand, with decreasing d, the i de- 
pendence of the susceptibility becomes much more com- 
plicated: The response below the second adatom (see 
the second maximum of |Xi U i b |) becomes stronger and 
stronger, the response at substrate sites between the 



adatoms increases and its absolute value develops a pro- 
nounced minimum close to 12, whereas the response be- 
yond the second adatom, for i > 12, gets very weak. 
Furthermore, while the susceptibility changes sign be- 
tween nearest neighbors, its two-site average is negative 
between the adatoms and also beyond the first one for 
i < it but is found to be positive for i > i 2 . There is 
another subtle observation, namely the (ferromagnetic) 
response at the nearest neighbor to the right of i\ is larger 
than the one to the left of i\ for all d down to d = 1, ex- 
cept for d = 5 and d = 9. The ratio x|"+x,i/Xi"-x,i is 
decreasing with decreasing d becomes smaller than unity 
for d = 5 and d = 9 and larger than unity again for d = 1. 

All these non-trivial features are perfectly captured by 
the R-DMFT and in fact result from an effective adatom- 
adatom interaction. This becomes obvious by comparing 
the results for the TIAM with those of a corresponding 



single-impurity Anderson model where the second (right) 
adatom a = 2 is missing or, equivalcntly where the hy- 
bridization to the second adatom is switched off. We have 
performed DMRG calculations for corresponding single- 
adatom models. To ensure a singlet ground state at half- 
filling, however, the substrate chain has to be shortened 
by one site on the right edge (L = 49). The resulting sub- 
strate susceptibilities xTi* are shown in Fig. [2] as green 
lines. 

Comparing the SIAM and the TIAM results to each 
other once more demonstrates that the effects of the in- 
direct nonlocal RKKY coupling become more and more 
pronounced with decreasing d. The differences between 
the single-adatom and the two-adatom physics visible in 
the susceptibilities for i > %2 are larger by more than an 
order of magnitude than the differences between the R- 
DMFT and the DMRG results. Again this shows that 
there are sizable effects on nonlocal magnetic correla- 
tions which do not fully feed back to the one-electron 
self-energy. 



B. Spin correlations and nonlocal susceptibitlities 

The DMRG data for x^" 13 and also for the equal-time 

spin-spin correlation function (S[S'[) are shown in Fig. [3] 
for d = 13 on a larger scale. Let us discuss the physics of 
this situation in detail. The response of the substrate to 
a static local field at (3 = 1 is governed by the low-energy 
excitations around the Fermi edge, i.e. u = 0. Contrary, 
the equal-time spin-spin correlation is obtained by a fre- 
quency integration of the dynamic (retarded) susceptibil- 
ity x™ b (u;) and thus includes several energy scales. Nev- 
ertheless, the spin-spin correlation behaves qualitatively 
very similar to Xi U i b j an d we will refer to this on an equal 
footing with the susceptibility. 

The ground state of the whole system is a spin singlet 
in all calculations discussed here. This provides us with 
a simple sum rule for the spin-spin correlation: Exploit- 
ing rotational symmetry, we have (S tot S{) = 3(Sl ot S{ z ) 

where S tot = S{ + S f 2 + S suh is the total spin and 
S sub = ^f =i fife the tQtal substrate spin . Usin g M tot = o 

in the ground state, we immediately find: 



(S{S{) + (S f 2 S{) + (S suh S{) = 
For the single-adatom model we then have 

(s{s{) 



{s suh s{) 



. 



(22) 



(23) 
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FIG. 3: (Color online) Spin-spin correlation function (SfSf) 
(upper panel) and magnetic susceptibility xT?i (lower panel) 
for U = 8 and V = \Pl as obtained by DMRG for a system 
with L = 50 (L — 49) substrate sites and two adatoms (one 
adatom) as functions of the substrate site i. Red lines: results 
for two adatoms and L — 50, TIAM. Green lines: results for 
L = 49 and a single adatom at the same position as the /3 = 1 
(left) adatom in the two-adatom model, SIAM. The dashed 
lines indicate the positions i\ and %2 of the substrate sites 
"below" the adatoms. 



but somewhat reduced as compared to the single-adatom 
model. 

Qualitatively the same applies to the susceptibility as 
can be seen from the lower panel in Fig. [3] For any large 
but finite system with a non-degenerate singlet ground 
state, we again have a simple sum rule: A singlet ground 
state and a finite gap implies that the total magnetic 
moment must vanish for any h\ up to some finite critical 



field: 



f\ 



sub\ 



0. Taking the 



derivative with respect to hi then yields: 



Xn +X21 



. 



(24) 



which explains why the spin-spin correlation is mainly 
negative: Namely, if summed over all substrate sites, 
it just compensates the adatom local moment. For the 
two-adatom model at d = 13, the adatom-adatom spin 
correlation (SlS^) is negative but its absolute value is 
small compared to (S{s[). Looking at Eq. (j22j) . the 
overall substrate response is thus still antiferromagnetic 



In the same way as above, ^ Xii h = ~Xii < for a 
single adatom, and for two adatoms the total response 
of the substrate is still negative but slightly reduced in 
absolute magnitude due to the presence of the second 
adatom since X21 < at d = 13. 

As can be seen in Fig. [3] by comparing with the re- 
sults for the single-adatom model, the most pronounced 
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effect due the presence of the second adatom consists in 
the strong enhancement of Xinii i- e - the response below 
the second adatom. This can easily be understood by 
referring to the RKKY limit for a system of finite size 
L: For V — > keeping U S> t fixed, charge fluctuations 
vanish and we are left with a Kondo-type model. In the 
weak-coupling limit J — > the substrate degrees of free- 
dom can be integrated out, and the adatom magnetic 
response, i.e. Xn an d X21, is perfectly described by an 
effective RKKY two-spin model 

-Hrkky = — JrkkyS{S2 > (25) 

where Jrkky = -f^Xi^ is given in terms of the static 
substrate susceptibility at J = 0. The Kondo effect, on 
the other hand, does not interfere with this picture as it 
is cut by the finite-size gap: One can define a coupling 
strength J c at which the Kondo temperature Tk becomes 
comparable with the finite-size gap. Then, for J < J c the 
Kondo effect is absent as there are simply no states at the 
Fermi energy available to screen the adatom moment^ 
This implies that the substrate is in a singlet state for 
weak J and thus x\ uh = J2iXa° = (d/dh^S? b ) = 0, 
i.e. there is no substrate contribution to the magnetic 
moment induced by the field at /3 = 1. From the sum 
rule Eq. ([2"4")l we thus have xn + X21 = 0, i.e. also the 
two adatom spins form a perfect singlet consistent with 
Eq. (|2~5j) . Hence, applying a field hi at adatom /3 = 1 
induces antifcrromagnctically aligned magnetic adatom 
moments with the same absolute magnitude. For J be- 
yond but close to the RKKY limit we therefore expect 
the absolute magnitude of the substrate response at i\ 
and «2 as almost equal. For finite and actually interme- 
diate J, see Fig. [3J the effect is strongly diminished but 
still clearly visible. Note that the above argumentation 
can analogously be given by referring to the spin-spin 
correlation. 

As mentioned before, looking at the sum rules ([221) 
and ((24)) . wc can understand that the response of the 
substrate is somewhat attenuated in the TIAM as com- 
pared to the SIAM. This reduction, however, is not ho- 
mogeneous: There is a comparatively strong reduction 
beyond the second adatom for i > ii while the response 
is nearly the same or even somewhat enhanced close to 
%2 for i < 12, and there is almost no effect for i < i\. 
That the effect is least pronounced close to the first 
adatom, can easily be understood by referring to the ex- 
treme Kondo limit where a picture of two separate Kondo 
clouds applies. In this case the magnetic response to the 
field applied to the first adatom would be the same as 
the response in the corresponding singlc-adatom model. 
As is seen in Fig. [3J however, close to i\ there are finite 
differences, i.e. the Kondo clouds do overlap, but the dif- 
ferences are small. Since according to the sum rule the 
total response must be weaker in the TIAM, a reduced 
response must and in fact does show up away from i\, 
i.e. for i > %2- 

The sum rule (|23)) for the SIAM may also be used 
to roughly estimate the size of the individual "Kondo 
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FIG. 4: (Color online) Susceptibility xfii a t the lattice site 
ii below the first adatom as a function of the adatom-adatom 
distance d for U = 8, V = y/2 and L = 50 substrate sites as 
obtained from R-DMFT (blue line) and DMRG (red line). 



clouds". Using the DMRG data for (S{S?), we define an 
integrated spin-spin correlation functional 

0(r)=l+ Y ^T"> (26) 
..ft- (S{S{) 

for the single-adatom model. We have (9(0) = 1. With 
increasing r more and more substrate spins around i\ are 
included in the sum, and 0(r) essentially decreases with r 
until — if all spins are included as is obvious from the 
sum rule Eq. (|23[) . The quantity gives the fraction of the 
adatom spin that remains unscreened by the substrate 
spins up to distance r from i\. Using a 90%-screening 
criterion, for example, i.e. 0(£k) = 0.1, the extent of 
the cloud amounts to £k — 10 — 15 lattice sites. This is 
consistent with the discussion given above. 

A criterion based on Eq. ([2^]) cannot precisely define 
the parameter range in which R-DMFT gives reliable re- 
sults. Fig.|4]demonstrates that, using R-DMFT, the devi- 
ation from the numerically exact DMRG data grows grad- 
ually when decreasing the distance between the adatoms 
d. 



C. Distance dependence 

To estimate the reliability of the mean-field approach, 
we focus on the susceptibility xli h &t the substrate site 
below the first adatom i = i\ where, according to the 
results shown in Fig. [21 the deviations are the strongest. 
Xffi is shown in Fig. [4] as a function of the distance d. 
There is a nice quantitative agreement of the R-DMFT 
with the exact DMRG result for large d. For smaller 
d, R-DMFT still predicts the correct trend, except for 
d= 1. 

For the same set of parameters Fig. [5] shows the local 
adatom susceptibility xu an d nonlocal adatom-adatom 
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FIG. 5: (Color online) Local adatom susceptibility xn and 
nonlocal adatom-adatom susceptibility \2l as functions of the 
distance d between the adatoms for U = 8, V = y/2 and 
for a system with L = 50 substrate sites as obtained by R- 
DMFT and DMRG (red lines). R-DMFT calculations are 
done with different numbers of bath orbitals in the effective 
single- impurity models: n s — 6, 8, 10, as indicated. 



susceptibility X21 as functions of the distance d = 4n + 1 
with integer n. In both cases the agreement of the R- 
DMFT with the DMRG results is excellent. Significant 
differences are found for d = 1 only and rapidly diminish 
with increasing d. 

Fig. [5] includes R-DMFT results obtained with differ- 
ent n s . On the scale of the figure, there is no difference 
between the results for n s — 8 and n a = 10 bath sites 
in the effective single-impurity model while the results 
obtained for \2\ with n s = 6 slightly deviate for inter- 
mediate distances around d = 25. This comparison shows 
that the R-DMFT results are converged with respect to 
n s . The differences to the DMRG data are thus intrinsic 
to the dynamical mean-field approach itself and not at 
all caused by discretization errors of the Lanczos solver. 

It is worth to mention that the distance dependence 
of X21 cannot be explained by conventional RKKY the- 
ory. For J — > 0, the magnetic susceptibility is deter- 
mined by the effective two-spin Heisenberg model Eq. 
(J2SJ which yields %2\ = ~Xu ~ 1/Jrkky with Jrkky oc 
(— l) d /d = 1/d at odd distances d. The decreasing ab- 
solute magnitude of X21 with increasing d and also the 
fact xn + X21 7^ just reminds us that with U = 8 and 
V 2 = 2 the system is well beyond the perturbative-in- J 
regime and that there is a strong substrate contribution 
J2iXa h necessary to fulfill the sum rule Eq. (|24|) . 

For large d the trends can rather be understood in 
a picture of two independent Kondo effects. Clearly, 
|X2i | is expected to decrease with d. More interesting is 
the behavior of Xn which develops a maximum around 
d = 15-20. The increase of xn with d at short d re- 
sults from a reminiscence to the RKKY limit: With in- 
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FIG. 6: (Color online) Local and nonlocal adatom suscep- 
tibilities xn and xn as functions of distance d for U = 8, 
V — a/2- R-DMFT calculations for different system sizes 
L as indicated. DMRG data are shown for comparison at 
nearest-neighbor distance d = 1 only. 

creasing d the effective coupling between the adatom de- 
creases and their magnetic moments tend to become free 
resulting in a higher local susceptibility xn- Substrate 
contributions are sizable but cannot outweigh this effect 
in xn, in contrast to X21 which is 2-3 times smaller in 
absolute magnitude. In the largc-d limit, where adatom- 
adatom interactions can be disregarded completely, one 
would naively expect a saturation of the local susceptibil- 
ity at the inverse Kondo temperature since xn °c 1/Tk 
in a single-impurity models However, xu must decrease 
since with increasing d at fixed L = 50 the adatoms move 
to the chain edges where we have a site-dependent Kondo 
temperature. This increases with decreasing distance to 
the edge as the non-interacting substrate local density of 
states at the Fermi energy is increasing. 

This interpretation is corroborated by Fig. [5] which dis- 
plays R-DMFT results for xn and X21 for different sys- 
tem sizes L = 30, L = 50 and L = 90. We find the same 
qualitative behavior in all three cases. Quantitatively, 
however, there are sizable differences at inter-adatom dis- 
tance d = 1, for example, which show that even with 
L = 90 substrate sites the chain center cannot be re- 
garded as bulk-like and that the center local density of 
states is still considerably dependent on L. On the other 
hand, the susceptibilities for d close to L, i.e. for systems 
with adatoms located at or very close to the chain edges, 
are almost converged. Note that xn for d = L is almost 
the same for L = 50 and L = 90. Again this shows that, 
at least for the larger systems, the magnetic response is 
dictated by the physics of the single-site Kondo effect, 
i.e. the presence of the second adatom has almost no ef- 
fect on xn and on the Kondo temperature of the first 
adatom. This does not exclude a finite magnetic interac- 
tion between the adatoms and in fact a non-zero X21 for 
d = L is found which, in addition, also does not depend 
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FIG. 7: (Color online) Local and nonlocal adatom suscep- 
tibilities xn and X21 as functions of V 2 /U for [7 = 8 and 
d = 1 as obtained by R-DMFT, real-space two-site DMFT 
and DMRG for a system with L = 50 substrate sites. Insets: 
same quantities plotted on a logarithmic scale. Dashed lines 
indicate the critical V 2 /U where xn and X21 diverge. 



on L for the larger systems. 

These adatom-adatom magnetic interactions are cor- 
rectly captured by the R-DMFT. R-DMFT and DMRG 
results coincide on the scale of the plot except for d = 1 
(DMRG results are shown for d = 1 only). This almost 
perfect agreement can be understood by referring to the 
strong-coupling limit: For J — > 00, the two adatom mag- 
netic moments form perfectly local Kondo singlets that 
do not interact with each other. This limit is trivially ac- 
cessible by the mean-field approach. For finite coupling, 
second-order perturbation theory in t/J predicts spin- 
spin correlations to decay as l/rf 2 j 68 i 69 The inter-adatom 
magnetic interaction is thus expected to scale as oc J~ 4 . 
This is also accessible to the R-DMFT approach while 
the neglected feedback of this effective interaction on lo- 
cal physical properties at one adatom, e.g. on Tk and 
thus on Xiii is °f higher order and small in the strong- J 
limit. 



D. Dependence on the local exchange coupling 

The breakdown of R-DMFT can be enforced, however, 
by decreasing J. Fig. [7J shows the susceptibilities for 
d — 1 as a function of V 2 /U. There is again excel- 
lent agreement for strong V 2 /U even with the simplified 
two-site R-DMFT. Deviations of the two-site approach 
from the exact xn and X21 become sizable for couplings 
smaller than V 2 /U « 0.5. The full R-DMFT is reli- 
able down to smaller values for V 2 /U but finally also 
starts to significantly deviate from the DMRG data for 
V 2 /U < 0.2. Here, as compared to the strong-coupling 
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FIG. 8: (Color online) Local and nonlocal adatom susceptibil- 
ities xn an d X21 as functions of V 2 /U for d = 1 and L — 50 
as obtained by R-DMFT. Results for different U as indicated. 
Upper inset: inverse susceptibility 1/xn and extrapolation 
(dotted lines) to 1/xn — 0. Lower inset: "phase diagram", 
separating the Fermi liquid state (FL) from an (artificial) an- 
tiferromagnetic state (AF) which shows up for weak V and 
strong U. 



limit, the local susceptibility is by more than an order of 
magnitude higher, i.e. the Kondo temperature is by more 
than an order of magnitude smaller (see the upper inset). 

For even smaller couplings, the mean-field approach 
breaks down completely and fails to maintain a Fermi- 
liquid ground state: The small- J limit is problematic for 
R-DMFT as the screening of the magnetic moments is too 
weak to compensate the ordering tendencies induced by a 
comparatively strong inter-adatom interaction. The sys- 
tem becomes too susceptible to an artificial spontaneous 
symmetry breaking that is induced by the mean-field ap- 
proximation itself. While the adatoms' state is given by 
a nonlocal SU(2) invariant singlet (| - \l^))/\/2 for 
J — > 0, the mean-field theory predicts an incoherent mix- 
ture of degenerate ordered states | fi) and 1 4-t) • 

This qualitative failure is indicated by divergencies of 
Xn and X21 which take place at coupling strengths V 2 /U 
that are somewhat smaller than those where first quan- 
titative deviations from the exact data were found (see 
insets in Fig. [7]). This also implies that the mean- field 
approach is able to exhibit its limitations by itself. 

Fig. [5] shows the susceptibilities as obtained by R- 
DMFT for d = 1 and different U and V as functions 
of V 2 /U. We find xn and X21 to diverge at the same 
point in parameter space. From extrapolations of the in- 
verse local susceptibility to 1/xn = at different U and 
V, shown in the upper inset, one may derive a mean-field 
"phase diagram" . This is displayed in the lower inset. A 
normal Fermi-liquid ground state found for large V and 
small U is separated from the SU(2)-symmctry-brokcn 
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antifcrromagnetic state realized for small V 2 /U. This 
"phase-transition line" should actually be interpreted as 
a crossover from the Kondo to the RKKY regime or, 
more precisely, as the boundary up to which R-DMFT is 
reliable. While the critical coupling V 2 /U is almost inde- 
pendent of U, there is some dependence on the distance 
d. However, this is weak: While a Fermi-liquid ground 
state is obtained down to V 2 /U « 0.14 for d = 1, we find 
a slightly smaller critical value of V 2 /U « 0.12 for d = 5 
and V 2 /U w 0.09 for d = 49 (at U = 8). 



E. Different distances between the adatoms 

More important for the reliability of R-DMFT is the lo- 
cal Kondo temperature. This becomes obvious if the two 
adatoms arc placed at a distance d = An + 3 with integer 
n, i.e. d = 3, 7, 11, ... etc. At the corresponding substrate 
sites i\ and ii (symmetric to the chain center) we have a 
low weight \Ui a k F \ 2 of the one-particle energy eigenstate 
of the non-interacting substrate at the Fermi wave vec- 
tor kp while |?7i Q fc F | 2 is high for distances d = An + 1. 
This pronounced odd-even effect is a consequence of sur- 
face Friedel oscillations. The weight |t r i Q fc F | 2 determines 
the local substrate density of states and thus also the lo- 
cal Kondo temperature. Consequently, Tk is small for 
d = An + 3 and the nonlocal RKKY interaction much 
more efficient. Using DMRG, we in fact find \x2i\ at 
d = 3 to be more than an order of magnitude larger 
than at d = 1. This effect even increases with increasing 
d — An + 3 since \Ui a k F \ 2 is decreasing if the adatoms 
move towards the chain edges. At the edges (d = 47) 
|^i a fep| 2 is suppressed by more than a factor 100 com- 
pared to the d = An + 1 case (d = 49), and the Kondo 
temperature is essentially vanishing. This regime is not 
accessible to R-DMFT. While the mean-field approach 
predicts the correct sign, the absolute value \x2i\ an d 
also Xn is strongly underestimated for d = An + 3. De- 
viations from the DMRG results grow with increasing d. 
This had to be expected, as a huge nonlocal suscepti- 
bility |x2i | induces via the Schwinger-Dyson equation a 
sizable contribution to the nonlocal self-energy which is 
neglected in R-DMFT. 

The ferromagnetic case is different. Here we consider 
distances d = 2n with integer n. In the Kondo limit of the 
model and for weak J, the ferromagnetic RKKY coupling 
of well-formed spins 1/2 leads to a triplet ground state 
as is easily verified by means of DMRG calculations for 
L = 50, i.e. there is a nonlocal spin 5 = 1. In all our 
calculations this spin is not screened by the substrate 
electrons. This may be explained by the fact that L = 
50 is still too small to accommodate the corresponding 
screening cloud. However, there is a nonlocal spin 5 = 1 
not only in the weak- J limit. In fact, we find a triplet 
ground state for any choice of U > and V ^ 0. 

Our R-DMFT calculations reproduce the spin-triplet 
ground state for small V 2 /U by predicting, for an in- 
finitesimally small external magnetic field in +z direc- 



tion, a spontaneously symmetry-broken ferromagnetic 
state | tt)- This corresponds to the M = 1 state of 
the DMRG spin triplet. With increasing V, however, the 
expectation value of the z-component of the total spin 
s tot = gf + s f + gsub deviates from unity andj beyond 

a critical hybridization V, even vanishes: (5* ot ) = 0. 
Hence, in the ferromagnetic case, R-DMFT is reliable 
in the small- J but appears to fail in the strong-coupling 
limit. This requires further investigations which, how- 
ever, are beyond the scope of the present paper. 



VI. CONCLUSION 

Conventional (RKKY) theory of indirect magnetic 
exchange predicts an effective exchange interaction 

^RKKYjj = J 2 Xif U ( w = 0) where Xij( u = °) is the non- 
local static susceptibility of the metallic host. This in- 
teraction survives, as a nearest-neighbor coupling^ even 
in the case of two magnetic impurities embedded in an 
infinite-dimensional lattice and is thus accessible by dy- 
namical mean-field theory. In the limit of infinite spatial 
dimensions or, at finite dimensions, within the dynam- 
ical mean-field approximation, one can therefore expect 
a finite response at one magnetic impurity subject to a 
local magnetic field at the other one, located at nearest- 
neighbor but also for larger distances. 

On the other hand, nonlocal effective interactions do 
not contribute to the single-particle self-energy on the 
DMFT level: The DMFT self-energy is just defined as the 
sum of the local skeleton diagrams only. This is a well- 
known shortcoming of mean-field theory which gives rise 
to artifacts in the RKKY limit. Namely, for J — > the 
magnetic impurities are only weakly coupled to the host 
and thus become extremely susceptible. A tiny Weiss 
field within DMFT is then sufficient to drive the system 
to an artificial symmetry-broken state, i.e. an antifcrro- 
magnetic state rather than a nonlocal singlet of the impu- 
rity magnetic moments is formed. A state with the char- 
acteristic distance dependence of the RKKY interaction, 
e.g. Jrkky oc 1/d for a one-dimensional system at half- 
filling, cannot be recovered within DMFT as it is always 
preempted by spontaneous symmetry breaking. Clearly, 
solutions with a finite magnetic moment could easily be 
suppressed in a mean-field approach. One should note, 
however, that the resulting magnetic susceptibility is un- 
physical, i.e. negative, as it refers to a thcrmodynamically 
unstable state. Therefore, in any case, the physics of the 
RKKY limit is not accessible by DMFT. 

The present study has shown, however, that beyond 
the RKKY limit, (real-space) DMFT is well suited to 
study even quantitatively the effects of indirect magnetic 
exchange. Here, we have concentrated on two magnetic 
"adatoms" on a one-dimensional "substrate surface" - a 
minimal model to study indirect magnetic interactions 
in competition with the Kondo effect for magnetic atoms 
on metallic surfaces and a model that is amenable to an 
exact numerical solution by means of the density-matrix 
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rcnormalization group. 

DMRG has been used to compute spin-spin correla- 
tion functions and static spin susceptibilities, and par- 
ticularly adatom-substrate site correlations and suscep- 
tibilities. Depending on the distance d between the two 
adatoms and depending on the hybridization strength 
V, rather complicated profiles are obtained. Comparing 
the results for the two-adatom (two-impurity) Anderson 
model with those obtained for the corresponding single- 
adatom (single-impurity) Anderson model, one can easily 
classify the different features of those profiles as single- 
impurity effects or as resulting from the adatom-adatom 
effective interaction. In this way, clear reminiscences of 
the RKKY interaction, i.e. of nonlocal singlet formation, 
are found to compete with the formation of Kondo clouds 
and screening of the adatom magnetic moments. In addi- 
tion, the profiles are strongly affected by the finite system 
size (chains with typically L — 50 have been considered 
here) and by effects resulting from strong surface Friedel 
oscillations in the local density of states, especially if the 
adatoms are in the vicinity to one of the chain edges. 
This complex interplay of different physical mechanisms 
is almost perfectly recovered by the real-space DMFT. 

Qualitatively, the real-space DMFT is reliable as long 
as the model parameters, in particular the local exchange 
coupling J cx V 2 /U, are in a regime well separated 
from the artificial symmetry-broken state. This param- 
eter regime, where the adatom susceptibilities are not 
too large or where the adatom magnetic moments are 
predominantly interacting with the substrate moments 
rather than among each other, however, goes well be- 
yond the extreme Kondo limit of non-overlapping Kondo 
clouds. The critical value for V 2 /U ~ 0.14 in units of 
t at d = 1 gives an impression of a lower bound for the 
applicability of R-DMFT. 

While the present study has focused on a one- 
dimensional model to allow for benchmarking against 
numerically exact DMRG results, future applications of 
the R-DMFT should address higher-dimensional systems. 
With increasing coordination number of correlated sites, 
the parameter space accessible to the mean-field ap- 
proach is expected to be become larger or mean-field ar- 
tifacts less pronounced. For the case of atoms trapped 
in optical lattices, there are impressing examples where 



R-DMFT has contributed to an understanding of the 
physics of inhomogencous systems with 0(100) corre- 
lated and geometrically inequivalent sites in two dimen- 
sions, for example^ 

Let us also point out that for Anderson-type multi- 
impurity or lattice models, one typically expects non- 
local magnetic correlations to diminish rapidly as the 
electron filling on the correlated sites is changed away 
from half-filling. Systems off half-filling are thus expected 
to be more amenable to an R-DMFT approach. At the 
same time they are also interesting physically as reduc- 
ing the filling away from half-filling affects local-moment 
formation as well. Hence, the competition between non- 
local RKKY interaction and Kondo screening must be 
seen as strongly filling dependent. 

For complex magnetic nanostructures with several 
magnetic adatoms in different chain or cluster geome- 
tries on two- or on semi-infinite three-dimensional metal- 
lic surfaces, a mean-field approach is inevitable anyway. 
Here the conceptual simplicity of a single-site mean-field 
theory, as compared to different possibilities for cluster 
extensions, is important as it allows to study almost ar- 
bitrary geometries. As in ab-initio studies, the accessible 
system size strongly depends on the remaining e.g. lateral 
spatial symmetries, and the computational effort scales 
nearly linearly with the number of inequivalent correlated 
sites only. 

The two-impurity one-dimensional Anderson model 
represents a model that is rather unfavorable to a single- 
site R-DMFT approach. Even for this the present 
study has shown, R-DMFT can in fact almost quantita- 
tively predict the effects of indirect magnetic exchange in 
competition with the Kondo and with geometrical effects 
- as long as the approximation predicts a Fermi-liquid 
ground state. 
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